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We examine an interacting dark matter-variable vacuum energy model for a spatially flat 
Friedmann-Roberston-Walker spacetime, focusing on the appearance of cosmological singularities 
such as big rip, big brake, big freeze, and big separation along with abrupt events {infinite 7 - 
singularity and new w-singularity) at late times. We introduce a phenomenological interaction 
which has a nonlinear dependence on the total energy density of the dark sector and its derivative, 
solve exactly the source equation for the model and find the energy density as function of the 
scale factor as well as the time dependence of the approximate scale factor in the neighborhood 
of the singularities. We describe the main characteristics of these singularities by exploring 
the type of interaction that makes them possible along with behavior of dark components near 
them. We apply the geometric Tipler and Krolak method for determining the fate of time¬ 
like geodesic curves around the singularities. We also explore the strength of them by analyzing 
the leading term in some geometric invariants such as the square Riemann scalar and the Ricci scalar. 
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I. INTRODUCTION 

In 1998 astrophysical observations coming from distant 
supernovae-stellar explosions led to the inevitable con¬ 
clusion that our Universe is actually speeding up rather 
than slowing down ^ due to the existence of a repul¬ 
sive agent known as dark energy. Such outstanding find¬ 
ing was promptly confirmed by additional cosmological 
observations based on the measurements of cosmic mi¬ 
crowave background anisotropies, baryon acoustic oscil¬ 
lations, and power spectrum of clustered matter [a [3]. 
Despite the plethora of observational evidences, in favor 
of the present acceleration of the Universe, accumulated 
over the last years by studying even more remote su¬ 
pernovae and launching new satellites [2], quite little is 
known about the true nature of dark energy. For in¬ 
stance, there is not a fundamental theory which can ac¬ 
count for the origin of dark energy at microscopic level. 

In order to shed some light over the dark side of the 
Universe one could devote some efforts to explore its ulti¬ 
mate fate [4]. An appealing route to follow is to analyze 
some dark energy models which can explain the current 
speeding up of the Universe but exhibit the presence of 
a cosmic singularity in the asymptotic (remote) future 
m, 0; so one may link the present state of the Universe 
with its drastic final state. A fundamental task to be 
addressed by these alternative scenarios is the classifica¬ 
tion of singularities that could emerge at a finite time or 
abrupt events [4], [5], [6], [7], [8], [9], [lO], [H], [12], [l3|. 
Several useful ways to characterize such final doomsdays 
are based on geometric methods which examine the ex¬ 
istence of causal geodesics that cannot be extended to 
arbitrary values of their proper time (geodesic incom¬ 
pleteness) [13] or the possibility to show that geodesic 


curves can be extended beyond cosmic singularity 
[ID- Another manner consists in taking into account 
the behavior of curvature invariants near the singularity 
so the strength of singularities can be determined using 
the necessary and sufficient conditions obtained by Tipler 
m and Krolak m- These conditions are considerably 
useful provided help to classify singularities as strong and 
weak types, giving some insights on the magnitude of the 
tidal forces experienced by a co-moving observer toward 
the singularity. In fact, the Tipler definition requires that 
any object has its volume crushed to zero as the singu¬ 
larity is approached m whereas the Krolak definition 
is weaker than the Tipler version and is related with in¬ 
quiries on the cosmic censorship conjecture [19] . 

A serious approach for understanding the nature of 
cosmological singularities requires to take into account 
the behavior of dynamical variables which enter into the 
field equations. Therefore the blow up of the energy den¬ 
sity and the divergent behavior of the pressure seem to 
be of interest along with the behavior of geometric quan¬ 
tities such as the scale factor, the Hubble function and its 
derivatives. In this way, one must focus on certain phys¬ 
ical properties which make these kinds of singularities 
fairly distinctive among them and how such traits can 
determine the ultimate fate of the universe 0> 0: 

0 , 0 , m, m, [ 13 ]. At this point, it would seem 
crucial to explore a viable cosmological scenario where 
the aforesaid singularities can appear naturally in order 
to explore its physical outcome. To do that, we are go¬ 
ing to present some interacting dark energy models [20] . 
m with the presence of singularities (or singular event) 
and abrupt events. A natural question that could arise 
is what kind of phenomenological interactions do imply 
the existence of a future singularity or abrupt event. We 
must stress that our approach for studying cosmic sin- 
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gularities is considerably different from previous works, 
thus, we focus in an interacting dark energy model with a 
nonlinear interaction term and then analyze the ultimate 
fate of the universe. As a consequence of this method, we 
will firstly study the existence of singularity along with 
its traits without imposing any particular type of scale 
factor. However, we will obtain the approximate scale 
factor or the exact one, in the cases where it is possi¬ 
ble, to illustrate the model. Since we do not postulate 
or give the form of scale factor from the very beginning, 
it is clear that only some singularities will emerge in our 
model. In this way, we will construct a posteriori a scale 
factor associated with a cosmic singularity for a given 
interaction. Such expression will be useful for studying 
then the behavior of the Hubble function, its derivatives 
in addition to the behavior of the energy density and 
pressure in terms of the cosmic time. 

Our goal is then to consider an interacting dark 
matter-variable vacuum energy (VVE) framework and 
explore the new details coming from the appearance of 
cosmic singularities in the remote future (distant past) or 
abrupt/singular events. Taking into account that these 
viable models exhibit a current accelerating phase to¬ 
gether with a future doomsday, we are going to describe 
the physical traits associated with them. Further, we an¬ 
alyze several kinds of singularities, make some comments 
on the type of interaction supporting them, in particular, 
we focus on the behavior of both dark component near 
the singularity. Finally, we examine these singularities 
with the help of the Tipler and Krolak criteria. 

II. INTERACTION FRAMEWORK 

To study the interacting dark sector model we consider 
the spatially flat Friedmann-Roberston-Walker (FRW) 
metric ds^ = —dt^ + a^{t)dx'‘dxi, where a{t) is the scale 
factor and dx'^dxi is the line element corresponding to hy¬ 
persurface of constant time. The universe is filled with 
two perfect fluids, one accommodates as a matter com¬ 
ponent while the the other represents a VVE substra¬ 
tum. Both fluids are described by linear equations of 
state, having energy densities pm, Px and pressures Pm, 
Px, respectively. The total energy density p and the con¬ 
servation equation p -\- 3H(p + p) =0 for this interacting 
two-fluid model can be written as 


P — Pm H" Px-i 

(1) 

II 

1 

3 

TS 

3 

1 

(2) 


where the dot stands for derivative with respect to the 
cosmic time ’ = d/dt being H = d/a, the prime is 
the derivative with respect to scale factor ' = djdri = 
d/iHdt = d/dln(a/ao)^ and oq is some value of ref¬ 
erence for the scale factor. We have assumed that 
both interacting components admit an equations of state 
Pi = ili ~ l)Pi with i = {m,x} such that the constant 


barotropic indexes jm and satisfy the next condition: 
0 < 7 a; < 7m- After having solved the algebraic linear 
system of equations Q-iH), we are able to get Pm and Px 
as functions of p and its p derivative p': 


IxP + p' 

Pm — 5 

(3) 

7m - lx 


+ 

II 

(4) 

'Im Tx 


To complete the model we introduce an exchange of en¬ 
ergy in the dark sector in terms of a factorized interaction 
3HQ{p^ p',p). Concerning this aim, we split the balance 
equation ([^ as follows: 

Pm TmPm Q•: (b) 

p'x + IxPx = Q- ( 6 ) 

From Eqs. (§-@ , we arrive at a second order differential 
equation for the energy density 

P" + (7m + lx)p' + ImlxP = Qilm “ 7x), (7) 

that we will call “source equation” henceforth. 

The uniqueness of the solutions corresponding 

to the algebraic system equations ([^(1^ allows us to ex¬ 
tract some interesting conclusions of the model. In the 
procedure outlined above one can find the energy density 
p by solving the source equation Q for a given inter¬ 
action term Q, subsequently, the energy densities of the 
matter and VVE components are reconstructed by means 
of Eqs. (|^-Q, pointing that this procedure does not rely 
on the specific cosmological equations that govern the dy¬ 
namic of an homogeneous isotropic flat universe. Later 
on, we will investigate a concrete example for a given in¬ 
teraction term and obtain several general features about 
the existence of initial and final singularities in the in¬ 
teracting dark sector model without knowing of the scale 
factor. 

As a result of this approach, we have been reducing 
the interacting framework to an effective one-fluid model 
with energy density p = Pm + Px and total pressure 
p{p, p') = Pm +Px = —p — p'■ Comparing the later equa¬ 
tion with the effective equation of state of the dark sector 
p = (7 — l)p , we obtain its effective conservation equa¬ 
tion 

p -h 7P = 0, (8) 

where the effective barotropic index reads 

ilmPm + ^xPx) 

7= -. (9) 

P 

In calculating the scale factor of the homogeneous and 
isotropic flat universe, we adopt the Einstein field equa¬ 
tions, so that the dynamic of the effective one-fluid model 
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will be governed by the corresponding Friedmann con¬ 
straint, 


= p. 


( 10 ) 


We will investigate a universe which transits from an 
initial matter-dominated phase into a final era dominated 
by an unknown component that will be identified with 
VVE, the former component is associated with an initial 
singularity and the latter one with a final singularity or 
a possible doomsday of the universe. In our model the 
VVE has an equation of state of the form = —px, so 
its barotropic index vanishes = 0 ), and therefore the 
source equation ([^ the matter energy density ([^, and 
the VVE density 13 turn to be given by 


p” + ImP' 
P' 


= lmQ, ( 11 ) 

Px = P + (12) 

Tm 


So far we have been speaking of an exchange of energy 
between matter and VVE fluids without specifying its 
form. We now propose that the interaction term in the 
dark sector can be defined via the equation 


Here we suppose that the input function / only depends 
on p, p', and p for simplicity. Replacing this interaction 
term Q into the source equation Q, it leads us to a 
nonlinear differential equation 


the pressure, the effective barotropic index, the acceler¬ 
ation term, the matter and VVE densities, all of them 
depend linearly with the input function / to describe the 
whole contribution coming from the energy transfer in 
the dark sector (13). An interesting point in regard with 
the role played by the interaction can be understood by 
analyzing a regime where / gives the largest contribu¬ 
tion, neglecting p and c. In such regime, the previous 
dynamical quantities can be easily found as follows 

P' « -7m/, P ~ 7m/, 7 ~ Imf/P, 


a « -Jmafl2, Pm « -px ~ -/■ ( 20 ) 

In order to further motivate our results, let us assume 
that function f{t) —>■ /(tg) when t —J, tg so that in the 
limit case /(tg) = ±oo, blowing up at the finite value 
of the cosmic time <g. Then, we have that the cos¬ 
mic time derivative of the Hubble variable H = p' 12 « 
—7m//2 —t Too, the pressure p « y^/ —t ±oo, the 
effective barotropic index y « Imf!P —t ±oo, the accel¬ 
eration a « —^rnO'f/2 —)■ Too and the matter and VVE 
densities pm ~ / —t ±oo, px ~ —/ —>■ Too, indicating 
that all these quantities diverge as t ^ tg. This promis¬ 
ing approach then opens the possibility of producing an 
initial or final singularity at tg within the framework of 
interacting dark sector, so it will be useful for examining 
the ultimate fate of a universe. In the next section, we 
are going to generate Q by selecting the input function 
/ and examine its physical outcome. 


III. INTERACTION FRAMEWORK PRODUCES 
INITIAL AND FINAL SINGULARITIES 


+ 7mP = -7m/'- 


(14) 


Taking into account the first integral of Eq. (14), it al 


lows us to write the equation of state of the mixture 
p = —p — p', the effective barotropic index y = —p'/P, 
the acceleration of the universe a along with the matter 
and VVE densities (fT3 as follows 


A. General properties 

In this section, we wish to investigate the initial and 
final singularities of the universe when the interaction in 
the dark sector is generated by the function /, written 
as a power-law of the energy density 


P' = -7m(cTpT/), 


(15) 


/ = ap-", 


( 21 ) 


P= ilm- 1)P + 7m(c T /), 


7 = 7n 


P P 


Pm=C + p+f, Px = -{c + f). 


(16) 


(17) 


” =-^(37m - 2 )p-^(cT/), (18) 

a 6 2 


where c is an integration constant. From Eqs. (15)-(19), 


we observe that the p derivative of the energy density, 


which yields an exchange of energy associated with the 
nonlinear interaction term (13) 


Q = nap " ^p' 


( 22 ) 


where a is a coupling constant and n is a non-vanishing 
real number. Then the first integral (15) becomes 


p' = -7m(cTpTap "■). 


(23) 


To simplify our formulation, we slightly redefine the in¬ 
teraction ( 22 ) to get a scenario where the integration con¬ 


stant c vanishes and does not contribute to the model. 


(19) y \/0 insert the first integral (23) with c = 0 into Eq. (22) 


and we obtain the final interaction 


Q = -na'yrnP 


— n—1 


(p + ap“ 


(24) 
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which can also be written as 


Q = -n 


ImP 


(25) 


From Eqs. (13 1 , ( pTj ) and ( p^ , we obtain the solutions of 
the source equation 0 in the case of c = 0. This kind of 
interaction is very interesting because the effective equa¬ 
tion of state of the effective one-fluid p = —p — p' for 
c = 0 turns to be that of a Chaplygin or anti-Chaplygin 
gas depending on the sign of the coupling constant a, as 
was noticed in Ref. [20]. Summing up, Eqs. (15)-(19l can 
be recast as 


7 = -7m(p + ap”"), 

(26) 

P = (im - l)p + a 7 mp“". 

(27) 

7 = 7m [1 + , 

(28) 

^ 1 /Q o^ —n 

= .(37m ^)P 0 P ^ 

a 6 2 

(29) 

a a 

Pm = P H ~ 1 Px = ~ ■ 

„n „n 

(30) 


We are in condition to go a step forward for integrating 
(261 and obtaining the energy density as a function of the 
sea e factor 


P = 


—a + b 




l/(n+l) 


(31) 


where cq and h are both arbitrary integration constants. 
A singularity can be achieved if any of the physical quan¬ 
tities, such as the energy density (31), the pressure (|27|, 


the barotropic index (28) or the acceleration (29), van¬ 


ishes or diverges at a finite time tg. One way to meet this 
condition is by choosing the integration constant b so that 
the square bracket in Eq. (31) vanishes as the scale factor 
reaches the finite value a^, namely b = a(as/ao)^'’'"*^""''^^. 


Finally, the energy density (31), its p derivative and the 


effective barotropic index are given by 


p= \ a 




37m (fl-1-1) 


l/{n+l) 


(32) 


p = -ajr, 


(7) 


37m(n-l-l) 


(33) 


7 = 


1 _ JL 


(S) 


37m(n-l-l) 


(34) 


To make contact with previous results we will examine 
the behavior of several important quantities explicitly. 


We start with the case a > 0. We achieve that the energy 
density p —>■ 0 for a —)■ with n > — 1 provided a < Ug 
or in the complementary scenario where a —>■ 00 and 
n < — 1. The energy density blows up (p —>■ 00 ) under 
two different conditions: i- a —>■ Og with a > Ug and 
n < — 1 or ii- a —>■ 0 and n > — 1. Another branch to 
study corresponds to a negative coupling constant, a < 0 . 
In this case, the energy density p —)■ 0 as a —)■ Og with 
a > Og and n > — 1. On the contrary, the energy density 
becomes divergent (p —>■ 00 ) in the a —> Og limit if two 
conditions simultaneously holds: a < Ug and n < — 1 . 
In the last case, , we also have that the energy densities 
tends to the constant which is i nde ed an exact solution 
p —>■ Pc = (—[ see Eq. (26)]. It can be obtained 
in the a —> 0 limit for n < — I or in the case a —> 00 for 
n > — 1 [see Eq. @]. This constant solution gives the 
exact de Sitter scale factor, 


= exp 


(_Q,)l/(»+l) 


-I 1/2 


(35) 


The p derivative p' = 2H is negative or positive definite 
and therefore the Hubble variable is a decreasing or an 
increasing function of the time, respectively. An inter¬ 
estingly fact is that the barotropic index always diverges 
(7 —> 00 ) in the o —> Og limit for any value of n. The 
remaining quantities, namely the pressure, the acceler¬ 
ation, and the dark matter and VVE densities can be 
obtained as functions of the scale factor by combining 
the equations (27), (29), (30) along with 


By considering the second term of the first integral (26) 
as the leading contribution [see Eqs. it is pos¬ 

sible to find the energy density after integrating the ap¬ 
proximate conservation equation (p « —•\/ 3 a; 7 mP“"“''^/^) 
to arrive at 


73 


ain 


2[v-l) 


At 


2 ( 1 - 1 ) 


(36) 


Now the approximate pressure (27) is given by 

-1 1^-2 


P « Q;7„ 


73 


a7„ 


2(12-1) 


At 


(37) 


while the approximate scale factors is obtained by inte¬ 
grating the Friedmann equation (10) as 


o « a. 


I - 


2(12-1) 


73 


a7„ 


2(12-1) 


iva'yra 

We introduce the main parameter 

2(n + l) 
2n + l ’ 


At 


(38) 


(39) 


which will be useful for describing the different kinds of 
singularities in the near future (see Fig. 1). We define 
At = t — tg with tg being a hnite cosmic time. Notice 









































5 



FIG. 1: Plot of the main parameter v in terms of n indicating 
the range where the different kinds of singularities or abrupt 
events are achieved. 


that the parameter v vanishes for n = — 1 or diverges 
at n = —1/2, hence, the expansions (36)-(38) are not 


well defined in those cases and we are going to deal with 
them separately. In fact, we will be able to solve the 
whole dynamics of the model and give the exact scale 
factor for both values of ly. 


Combining the energy density (36), (24), and (26) 


(30), we can rewrite the interaction term, the approx¬ 


imate Hubble variable, its first, its second time deriva¬ 
tives, the acceleration term, the barotropic index, and the 
matter and VVE dens itie s, under the assumption that 
the ap“”-term in Eq. (26) becomes dominant: 




2{iy-2){iy-l) 

37,nAt2 


(40) 


H «-= 


1 


V3 


a7„ 


2 ( 1 .- 1 ) 


At 


I/-1 


(41) 


« - 


a7„ 


73 


Q:7„ 


2 (z.-l) 


At 


u-2 


(42) 


H 


73a7m(^- 2) 

4(1.- 1 ) 


73 


a7„ 


2 ( 1 .- 1 ) 


At 


u—S 


(43) 


- 


( 44 ) 


7 


aim 


2{u-l) 



(45) 


_ -2H _ 2H 

Pm ~ 5 Px ~ ■ 


(46) 


We have carried out a detailed calculation of the dynam¬ 
ical quantities (densities, pressures, interaction term, ac¬ 
celeration, etc) that we will employ for addressing the 
issue of cosmic singularities. The next step is to classify 
the singularities using the behavior of the above quanti¬ 
ties near them. 


B. Classification of singularities 


After having discussed the main traits of an interact¬ 
ing cosmology with two components in the presence of an 
initial or a final singularity, we must classify the different 
kinds of singularities that will emerge within this inter¬ 
acting framework. In order to do that, we mainly use the 
approximate energy density (36), the pressure (37), the 


scale factor (38) and remaining quantities (40)-(46). We 


classify the singularities according to the values taken by 
the main parameter, v. Here, we will exclude the non 
interacting cases corresponding to ;. = 2 and v = 1. 


1. {1 < V < 2)-case Big Brake singularity 

If we now consider the case of a positive coupling constant 
{a > 0) and the time approaching to the finite value tg 
from the left t ^ tg (thus At < 0) then the scale factor of 
the universe (38) reaches a finite value a —>■ Og, the time 
derivative of the scale factor vanishes a —)■ 0 but the sec¬ 
ond and third derivatives both diverge at the singularity 
(ds = —oo, dg = -|-oo). However, the energy density (36) 
along with the Hubble variable (41) are zero near the 


singularity, namely p{ts) = 0 and H{tg) = 0. The time 
derivative of the Hubble variable and their subsequent 
ones diverge as well as the acceleration of the universe 
d « UgH —)• — 00 , see Eq. (42), while the pressure (37) 


positively grows without limit {p —)■ -l-oo). Consequently, 
the interacting dark sector model has a late-time singu¬ 
larity at the finite cosmic time tg for n > 0 and a > 0 
characterized by a finite scale factor, a vanishing time 
derivative, energy density and Hubble variable as well. 
But the acceleration and the pressure both diverge while 
the effective fluid fulfills anti-chaplygin gas equation of 
state (27). These results show us that this kind of be¬ 


havior corresponds to a big Brake singularity, see [7] . In¬ 
terestingly enough, it turned out the energy transfer from 
the VVE to the matter diverges as t ^ tg (Q ——oo). 
Pm —>■ + 00 , and p^ —>■ —oo [see Eq. (40)]. It is important 
to stress that the physical set up associated with a big 
brake singularity can be described in terms of a tachyonic 
scalar field [7]. A weaker extension of the big brake sin¬ 
gularity is obtained when a dust component is included 
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in the Friedmann equation provided the Hubble function 
does not vanish any longer at the singularity [ 8 ]. 

We end this case with the analysis of a complemen¬ 
tary branch where the scale factor is defined in the re¬ 
gion t > ts and the coupling constant is negative (a < 0 ). 
The former branch which corresponds to a > 0 and t < tg 
matches with the latter one at t = tg up to its first time 
derivative whereas their higher order derivatives are all 
divergent [S]. Near the singularity the effects introduced 
by the interaction (40) on the matter and VVE compo- 
are the same of the a > 0 case, thus, their 


nents ( |30[ 
limits are coincident. 

2. {2 < ly < c»)-case Big Separation singularity 

In the interval of v € (2,oo) under the assumption of 
positive coupling constant (a > 0 ), we find that the scale 
factor (38) reaches a constant value a., as t —>■ tg from 


the left while the energy density (36) and the pressure 
(37) both approach to zero [5] , however, the barotropic 
index (45) becomes divergent as t —>■ tg. The first time 


derivative of the scale factor a —)■ 0 , its second derivative 
d —>■ 0, the Hubble variable H —>■ 0 and its first time 
derivative —>■ 0. Taking into account that p' = 2H —>■ 0 
ast ^ tg and p" = 2H/2>H, the source equation (11) then 
becomes 


H 


^HQ. 


(47) 


This equation explicitly shows that the second and subse¬ 
quent time derivatives of the Hubble variable are strongly 
dep enden t on the interaction term. In fact, from Eqs. 
(43) and (47), we find that H, H, and their subsequently 
time derivatives vanish up to order k < 1 ]] while 

the remaining ones diverge for k > [[ly — 1 ]] if t —>■ <s, 
where [[a;]] stands for the integer part of its argument. 
From Eqs. (40) and (46), we also have that Q —> oo. 
Pm 0 and px 0 as t —> ts- In conclusion, the 

interacting dark sector model exhibits a big separation 
singularity, sometimes considered as a more softer singu¬ 
larity. 


3. 


-|-oo)-case Exact solution 


We will discuss in detail the mysterious case associated 
with n —>■ —1/2. In fact, we are going to show that the 
field equations can be solved exactly and therefore we can 
obtain a complete picture of the interacting model. To 
this end, we introduce a new variable s = so 

that the Friedmann equation (10) with the energy density 


(48) 


(32) reduces to the linear differential equation 

s = ± - s). 

Having solved this equation we obtain four different types 
of solutions which need to be examined. Our starting 
point is to consider the expanding ai-solution and its 
effective barotropic index (7 = —2H/3H^), 


1 + 


At 


(49) 


7i = -7me 


— Ot'yrn‘^t 


(50) 


where i7i = cti/oi > 0 and 71 < 0. The contracting 03- 
solution, obtained from the time reversal symmetry of 
the former one, leads to 03 = ai{—t) with = 03/03 < 
0 and 73 = 7 i(—t). We name the another expanding 
solution as 04 and is given by 


04 = Os 1 — e 


-Oi-y-n 


74 = 


(51) 


(52) 


Here 774 = 04/04 > 0 and the contracting case is defined 
as 02 = 04(—t) with H2 = a2/(i2 < 0 and 72 = 74(—t). 
There we have assumed a positive coupling constant, 
a > 0 . The branches oi and 03 are defined in the in¬ 
terval —00 < t < 00 while for the branches 02 and 04 are 
specified in the intervals t < 0 and t > 0, respectively. 
We have fixed the integration constants so that the so¬ 
lution oi —)■ Os and 03 —)■ Os in the limits t —)■ —00 and 
t —> -|-oo, respectively. We have fixed the final big crunch 
or the initial big bang singularity in 02 or in 04 at t = 0. 

The solution oi (03) describes a universe which ex¬ 
pands (contracts) from a finite scale factor Os, free of an 
initial singularity, (from a infinite scale factor) at t = —00 
and a final de Sitter stage (a finite scale factor). However, 
the solution 04 (02) represents a universe which evolves 
from the big bang singularity (a finite scale factor Os) at 
t = 0 (in the limit t = —00) and increases (decreases) 
monotonically until reaches the finite scale factor o^ (the 
big crunch singularity) in the limit t —>■ 00 (at t = 0 ). The 
last two scale factors 04 and 02 behave as a power-law 
solution o « as[-\/3a7m |t|/2]^/^^"* in the limit |t| —0 
when the universe begins from a big bang singularity or 
ends in a big crunch singularity. 

The scale factors 03 and 04 (oi and 02) all of them go to 
ttg as At —7- 00 (At —>■ —c») while its first and higher order 
time derivatives vanish in those limits. Furthermore, we 
find that the Hubble variables can be written in terms of 
the barotropic index as 


H,, = 


1 




(53) 


where the index i runs from 1 to 4. Their subsequent 

(k) 

time derivatives 77/ ' —0 as a —>■ Ug. 

The energy density of the mix, its pressure and the 
interaction term read 


Pi — 


-l+(^^ 
a. 


Pi = (7m - l)Pi + O-lmpJ’^ ■, 


(54) 


(55) 


Qi = 


a7„ 


a + 7i 


1/2 


«! = a. 


(56) 
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At this point, we highlight that the time dependence of 
these quantities can be obtained by replacing each one 
of the exact scale factors found before [see Eqs. (49) 
and ([4^]. From (54) and (55), we have that the energy 
density p{a) —>■ 0 and the pressure p —?► 0 in the limit a —>■ 
Os. Although the energy density and the pressure vanish 
in this limit, however, the effective barotropic index 7 = 
ip + p)/p = 7 m[l + OLp~^l'^\ diverges. Most importantly, 
the higher order time derivatives of 7 also diverge as a —> 
tts so we name this behavior as “infinite 7 -singularity”. 

Consequently, the universe presents an infinite 7 - 
singularity in the distant past or the remote future 
(t = ± 00 ) where the scale factor reaches a finite val¬ 
ues Us, but the first and subsequent time derivatives of 
the scale factor vanish, further, the Hubble variable along 
with their higher order time derivatives also vanish near 
ttg. We remark this “new singularity” differs from the 
one reported in [53] provided it could not be reached 
at finite time. Nevertheless, this singularity is charac¬ 
terized by a non-vanishing interaction term (56) at Ug, 
Q —> a^ 7 m/ 2 , while the matter component pm 0 and 
the VEE px —>■ 0. Even though all physical quantities (p, 
p. Pm, Px ) fades away quickly as one approaches to the 
singularity, in the far remote future or distant asymp¬ 
totic past, the interaction cannot be turned off neither 
the barotropic index. 


4- (—00 <v < 0)-case New w-singularity 

In addition to the spacetime singularities at finite time 
mentioned above there are another kinds of abrupt events 
[ 22 ] related with the ultimate fate of the universe which 
indeed requires some mentioning for sake of complete¬ 
ness. The abrupt events were discovered within the 
framework of dark energy model and attracted some at¬ 
tention mostly beacuse are less dangerous than singular¬ 
ity. Concerning this aim, we are going to analyze a new 
kind of abrupt event that is related with the so called 
w-singularity found in |23j and re-examined in [24] . 

For these values of the main parameter v, the scale 
factor (38), the energy density (36) and the pressure (37) 
Og, p —>■ 0 and p —>■ 0 as the time 


For this particular value of the main parameter v, the 


have the limits a 
variable t —> 00 while the first and higher order time 
derivatives of these three quantities vanish in the remote 
futureft —» 00 ). Due to that the effective barotropic in¬ 
dex (45) and their subsequent time derivative 7 ^^^ di¬ 


verges up to order k < [[—but vanishes for all the 
following ones in the limit t —>■ 00 ; as a result the uni¬ 
verse ends in an abrupt 7 -singularity in the remote fu¬ 
ture, again, this kind of singular event differs from the 
case explored in [23] . The approximate interaction term 


(40) vanishes in the abrupt event whereas Eq. (46) shows 
us that Px ^ 0 and pm —>■ 0 as a —>■ Og. This singular 
event differs from the “infinite 7 -singularity” because the 
interaction vanishes in the latter case while the former 
one leads to Q —>■ Q;^ 7 m/ 2 . 

5. {v = 0)-case Big Rip/Crunch Singularity 


constant n = — 1 and the interaction term ( 22 ) becomes 


a linear function of the 77 -derivative of the energy density, 
Q = —ap', and therefore the first integral (23) turns to 
be 


p' + 7m(l + a)p = c'ym- 


(57) 


We will solve this equation for the energy density and 
below show the general exact solution of the Friedmann 
equation for any value of the integration constant c. This 
includes the corresponding power law scale factor for c = 
0 : 


= a,t2/37™(l+a) 


c = 0 , 


P = 


(?) 


37m(H-a) 


c = 0 , 


37m (1 + a)^ ^ « 

and two families of solutions associated with c ^ 0, 


1 — cos Lot 
2 


l/37„(H-a) 


P = 




UJ^ 


sin cat 
1 — cos UJt 


where -|- a) > 0 and 


a = 


1 -I- coshwt 


-I l/37„(l+a) 


P = 


3cW; 


.2^,2 r 


sinh Lot 
1 -|- cos uot 


(58) 


(59) 


(60) 


(61) 


(62) 


(63) 


where = —3c 7^(1 -I- a) > 0 while the pressure and 
the interaction term are 


P = C7m + [7m(1 + a) - l]p, 


Q = a7m[c-|- (1 -I- a)p]. 


(64) 


(65) 


For 


= 3 c 7^(1 -|- a) > 0 and a > — 1 , we have 


that the integration constant c > 0 and the solution (60) 


represents a universe with a finite time span that begins 
with a big bang at t = 0 , then passes by a maximum 
at tmax = tt/w and ends in a big crunch at tbc = 27r/a; 
[B]. However, for a < — 1 and c < 0 there is a significant 
difference in the behavior of the solution (60) provided 


the universe has a finite time span but begins with a 
contracting phase at t = 0 associated with an infinite 
scale factor, then bounces at tbounce = and ends 
with an infinite scale factor in a big rip singularity at 
tbr = 27r/a; [ 6 ]. Near the big bang, big crunch and big 
rip singularities, the energy density (61), the pressure 


(64), interaction term (65) and the matter and the VVE 
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densities (46 1 diverge, namely p —>■ oo, p —>■ —oo, Q 
+ 00 , Pm —t — oo, and px —> +oo. At the extrema of 
the scale factor (601 tmax/bounce = we have that 

p vanishes, both p and Q become constants p = cjm, 
Q = ap while Pm = Px = 0. 

For = — 3 c 7^(1 + a) > 0 and a > —1, we have 
that the integration constant c < 0 and the solution (62) 


In same manner, a spacetime is Krolak strong |18) iff as 
the proper time t ^ the integral 


describes a universe with an initial contracting de Sitter 
phase in the remote past, then bounces at t = 0 and ends 
in an expanding de Sitter stage in the remote future. For 
a < — 1 and c > 0 , the universe evolves from a zero radius 
(vanishing scale factor) at the distant past, then passes 
by a maximum at t = 0 and ends in the remote future, 
again, with a vanishing scale factor. The energy density 
([63|), pressure (1^, the interaction term (1^ and the 


matter and VVE densities (30) approach to the constants 
p p = -p,Q ^ a7m(l - c). 

Pm —t —c and Px —>■ ca/(l + a) as time t —>■ ±oo. At 
t = 0 the scale factor exhibits an extreme where p = 0, 
P = cjm, Q = ap and pm = Px = 0. 


6. {0 < ly < l)-case Big Freeze 

The big freeze singularity was firstly encountered in the 
literature by Bouhmadi-Lopez et al. within the context 
of a phantom generalized Chaplygin gas m- Further, 
they proved that this singularity at a finite scale factor 
arises in a Randall- S undrum I brane-world scenario if the 
brane is filled with the dual version of the generalized 
phantom Chaplygin gas HU. Besides, the avoidance of 
this singularity within the context of quantum cosmology 
was explored with the help of the Wheeler-de Witt equa¬ 
tion by mimicking a (phantom) generalized Chaplygin 
gas with a scalar field m- 

The scale factor (38) approaches to a finite value a —>■ 
Qs in the limit t ^ tg. However, its time derivative a and 
their subsequent ones with A: > 1 diverge for t ^ tg, 


indicating the existence of the big freeze singularity m- 
Concerning the Hubble variable (41), the energy density 


(36), the pressure (37), the barotropic index (45) and the 


interaction term (401, these and their higher order time 
derivatives diverge in the limit t —>■ tg. Then, we have a 
divergent behavior of the dark matter and VVE density 
as t —>■ A. 


C. Krolak and Tipler criteria 


In what follows, we give a reinterpretation of the above 
classification making a description of the singularities 
and/or abrupt events from a geometric point of view 
based on a method developed by Tipler m and Krolak 
[18]. A spacetime is Tipler strong m iff as the proper 
time t —>■ tg, the integral 

T{t) = f dt' [ \TZabU°‘v!’\dt" oo. ( 66 ) 
Jo Jo 


IC{t) = f \TZabU°'u^\dt —>■ oo, 
Jo 


(67) 


where the component of Ricci tensor are understood to be 
written in a parallel transported frame along the geodesic 
curves. Notice that a singularity can be strong by Krolak 
criteria but weak according to Tipler’s criteria, however, 
the reverse situation always holds. Because weak singu¬ 
larities can be extended beyond them, the method de¬ 
veloped by Tipler and Krolak are useful tools for deter¬ 
mining the fate of the universe in terms of the fate of 
geodesic curves near potential strong singular point. 

Let us consider time-like geodesic curves, x® = c 
with i spatial index and c a constant [53] ? associated 
with co-moving observer, i.e, we take into account a co¬ 
moving world-line congruence with velocity ri“ = (9t)“ = 
( 1 , 0 , 0 , 0 ) so that the proper time and the coordinate time 
are the same. Moreover, the components of the Ricci 
tensor measured by an observer along this congruence 
lead to TZ abU°'u^ = —3d{t)/a{t). Using the latter fact 
together (38) and (42) is not difficult to see that a big 


brake singularity is T-weak and K-weak provided both in¬ 
tegrals (66)-(67) vanish for i' G ( 1 , 2 ), namely T oc At'^ 
and JC oTAt’^ as At -g 0. Nevertheless, the leading 
term in the square Riemann 'R-°‘^'^'^'R-abcd cx (d/a)^, given 
by diverges as At —)■ 0 (23 . Also the Ricci 

scalar, TZ oc d/a oc At''~‘^, blows up at the singularity 
pSj . A big separation singularity can be considered as 
a weaker event provided is not only T-weak and K-weak 
but also TZ°‘^‘^‘^TZabcd and TZ both vanish at the singular¬ 
ity. The “new w-singularity” is T/K-weak and avoids 
divergences in geometric invariant as TZ'^^'^'^TZabcd in the 
asymptotic past or remote future, indeed a similar situ¬ 
ation occurs for the infinite 7 -singularity. On the other 
hand, the behavior of scale factor near a big rip singu¬ 
larity is a{t) oc (t — tbr)”^ with p > 0 , so we classify it as 
a K/T-strong singularity, and all scalar invariants blow 
up at tbr- Besides, one could expect that both big bang 
and big crunch singularity exhibit a similar behavior re¬ 
garding the K/T criteria or the blow up of TZ°’^‘^^TZabcd- 
In the case of a big freeze ultimate fate, the leading term 
in Krolak strength gives K, oc At'^~^ so it diverges as 
t ^ tg, however, Tipler measure involves a second inte¬ 
gration and yields T oc {t — tsY, being totally regular 
as t —>■ ts provided v G (0,1). Notice that the squared 
Riemann and the Ricci scalar both diverge for the big 
freeze event, as can be seen from (38) and (42). It is im¬ 


portant to stress that the analysis performed with casual 
geodesic which meets a T /AT-strong singularity gave the 
same kind of finding because involves the integral of com¬ 


ponent Riemann tensor parallel transported, TZ\^jvJu'^, 
for which the non-vanishing components turned to be the 
same {TZ’^tct = —d/a) also [15] ■ We end this section by 
studying the violation or not of the energy conditions 
for the aforesaid singularities [2U- Interestingly enough, 
we obtained that all the energy conditions (WEC, NEC, 
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SEC) are fulfilled near the singularities except for a few 
cases, namely the big rip and the infinite 7 singularity 
(associated with the expanding solution oi) both violate 
SEC [H]. 

An open question related with the present work is what 
kind of fundamental theory can support these singulari¬ 
ties? In this direction, Barrow and Graham have recently 
shown that a new ultra-weak generalized sudden singu¬ 
larity can be supported by a scalar field with a simple 
power-law potential [28] . We hope to return to this issue 
in the near future. 


IV. CONCLUSIONS 

We have investigated the evolution of an interact¬ 
ing dark sector model in a spatially flat FRW space- 
time, solved exactly the source equation when the en¬ 
ergy exchange is produced by a nonlinear interaction 
term and found that the equation of state of the effec¬ 
tive dark fluid turned to be that of a Chaplygin or anti- 
Chaplygin gas according to the sign of the coupling con¬ 
stant. We highlighted that this interacting dark matter- 
VVE model, generated by a nonlinear interaction term, 
harbored future-like singularities related with the ulti¬ 
mate fate of the universe and two abrupt events called 
infinite 7 singularity and new w-singularity. 

We have obtained the approximate scale factor in the 
vicinity of the singular events and used it to study physi¬ 
cally viable models describing the current along with final 


state of the universe, in particular, we have examined the 
behavior of the interaction and the dark component en¬ 
ergy densities near the singularities. We have classified 
entirely the singularities of the model for any value of 
the main parameter v while for certain value of we 
have shown that the two-fluid system is fully integrable, 
so the exact form of the scale factor and the remaining 
representative quantities of the model were found. 

Using the aforesaid scale factor and the geometric 
method developed by Tipler and Krolak for the case of 
time-like geodesic curves, associated with co-moving ob¬ 
server, we have obtained a regular behavior of Tipler 
measure in the case of big separation, big brake, big 
freeze, new w-singularity, and the infinite 7 -singularity, 
making them traversable in the sense that a time-like 
geodesic curves can be extended beyond such singular 
events. However, big freeze singularities are K-strong due 
to Krolak measure becomes divergent near them. On the 
other side, we have shown that brig rip singularity cannot 
be smoothed within interacting scenarios and therefore 
Tipler and Krolak measures diverges as geodesics meet 
such final event. 
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